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In grand unified theories (GUT) based on SO(10) all fermions of one generation are embedded 
^ : in a single representation. As a result, the top quark, the bottom quark, and the r lepton have 
a universal Yukawa coupling at the GUT scale. This implies a very large ratio of Higgs vacuum 
expectation values, tan/5 ~ mt/rrn,. We analyze the naturalness of such a scenario quantitatively 
including all the relevant radiative corrections and find that in minimal unified supergravity models 
with universal soft supersymmetry breaking parameters at the GUT scale, the necessary amount of 
fine-tuning needed is excessive. GUT threshold correction to the universal Higgs mass parameter 
can significantly reduce the fine-tuning required for such large values of tan (3. We also point out 
that the top quark mass can play a crucial role in explaining the hierarchy between the SUSY 
breaking scale and the electro-weak scale and, hence, the naturalness of large values of tan f3. 
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1. Motivation 

One of the main motivations for supersymmetry (SUSY) is the cancellation of quadratic 
divergences in unrenormalized Green functions. This stabilizes any mass scale under radiative 
corrections and allows the existence of the large hierarchy between the Planck-scale and the 
electro-weak scale. Furthermore, it has been shown recently that in the minimal supersymmetric 
extension of the standard model (MSSM) [1] the SU(3) C <g> SU(2) L <g>U(l) Y gauge couplings unify 
at a scale Mqut — 2 x 10 16 GeV [2]. Additionally, the unification of r and bottom Yukawa 
couplings at Mqut can be achieved within the MSSM for particular values in the m^-tan (3 plane 
[3]. Further unification can be achieved if we embed the SM gauge group within an SO(10) [4]. 
Here all the fermions of one generation fit in one single 16 dimensional spinor representation F(16)j 
(i = 1, 2, 3 is the generation index). Furthermore, SO(10) is anomaly-free and thus provides a nice 
explanation for the cancellation of anomalies within one generation of standard model quarks and 
leptons. 

If we assume that the two Higgs doublets, h u and hd, required in supersymmetric theories to 
give masses to up- and down-type fermions, are embedded in a single representation, h(10), and 
that all the Yukawa couplings of the third generation come from the unique coupling F3 £g> F3 ® h 
we find 

yt(M GlJT ) = y b (M GlJT ) = y T (M GVT ) . (1) 

In these models, the large ratio of mt/m^ can be explained by a large ratio of Higgs vacuum 
expectation values (VEVs) 

tan/3 ^ ^4^^ = 0(50). (2) 

The prediction for the bottom quark mass from r-bottom Yukawa unification tends to be larger 
than its experimental value, m& = 4.7±0.2 GeV [5], unless a significant part of the QCD correction 
is cancelled by a large top Yukawa coupling [3]. This predicts that the top quark mass, mt, must 
be close to its infra-red (IR) fixed point which is roughly obtained by setting the /^-function to 
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zero, i.e. 



m t m b _ /onn n^n2 



, „ ~ (200 GeVV . (3) 
sin 2 /? 6cos 2 P K J K ' 

For a top quark mass of m t = 176 ± 8(stat) ± 10(sys) GeV suggested by recent CDF results [6] 

eq. (3) has two ranges of solutions corresponding to tan/3 ~ 1 and to very large values of tan (3 

that are remarkably compatible with eq. (2) [7]. 

In minimal unified supergravity (SUGRA) models, where all the soft SUSY breaking mass 
parameters are universal, the Higgs parameters m\ (h — h u , hj) are driven to negative values by 
the large Yukawa couplings [8]. As a result, the correct pattern of electroweak symmetry breaking 
occurs naturally in a large region of the parameter space. The resulting ratio of Higgs vacuum 
expectation values can be expressed in terms of the parameters of the Higgs potential evaluated 
at the electroweak scale set by the Z boson mass, m z = 91.187 [5]. At tree-level the minimum 
conditions are 

tan2 , = <±4±l4 . (4) 

m h u + 2 m z 

sin 2(3 = - - . (5) 



m 



hd 



+ ml + 2^ 



where B denotes the running B parameter evaluated at the electro- weak scale. In the large tan (3 
limit we find that [i has to be tuned such that the denominator of eq. (4) vanishes as sin 2 2(3. 
The absence of a hierarchy problem then just means that fi 2 and m| should be of the order of 
or not much larger than m 2 . This feature is common to all SUSY models. The second condition 
[eq. (5)] only becomes problematic in SO(10) type models with large values of tan/5 [9]. It has 
been pointed out in ref. 10 that an approximate symmetry is the only case where the small values 
of B\i required in SO(10) by eq. (5) are stable under radiative corrections and thus natural. In 
this paper, we will combine such a symmetry with the radiative electro-weak symmetry breaking 
scenario in minimal SUGRA models and try to analyse the required fine-tuning quantitatively. 

Our paper is organized as follows: in chapter 2 we review the fine-tuning problem with uni- 
versal soft SUSY breaking masses. In chapter 3 we allow for non-universal soft SUSY breaking 



parameters at Mqut m order to accommodate an approximate symmetry that yields large values 
for tan j3 and in chapter 4 we will present our conclusions. 

2. The fine-tuning problem with universal mass parameters 

In minimal N — 1 SUGRA models the effects of SUSY breaking can be parameterized by only 
5 (4 in the MSSM) parameters 

V S oft = (AW 3 + BW 2 + CWi + H.c.) + m 2 ^ cfifa + m 1/2 ^ (^a + H.c.) , (6) 

i a 

where Wi, W2 and W3 are the linear, quadratic and trilinear parts of the superpotential, respec- 
tively (in the MSSM, W\ = 0). Furthermore, itiq (m^ 2 ) is the universal mass parameter of all 
the scalar fields 4>i (gaugino fields ip a ) at the Planck scale, Mp. In general, the universality of 
all the soft SUSY breaking parameters is violated through radiative corrections and significant 
non-universal terms can already be generated by renormalization group (RG) evolution from Mp 
to Mqxjt [11]- 111 this chapter, we will follow the common practice of neglecting these model- 
dependent effects. 

We will now investigate the feasibility of incorporating a symmetry that sets fiB = 0. Possible 
candidates for such a symmetry advocated in ref. 10 are the Peccei-Quinn (PQ) symmetry which 
implies that /i — and R symmetry [12] which implies that B = 0. The latter requires also that 
A = mi/2 = as can be seem by inspecting the f3 function of B. The amount of fine-tuning 
required to obtain a particular value of /3, can be defined as [13;14] 

B d sin 2/? _ B_dB_ _ B 
B ~sin2/5 dB ~BdB~B' 

For a natural solution we expect Rb = 0(1) and an increasing Rb means increasing fine-tuning. 
Since the set of RG equations for the soft SUSY breaking mass parameters [15] is homogeneous 
we can write the solutions as 

= C iI, 



(8) 

Here, the linear input and output parameters are denoted by I = m^^A^B and O — rrii (i — 



O 2 = C 2 OI I 2 . 



1,2,3), Af (f — t,b,r), B, respectively. The quadratic input and output parameters are denoted 
by I = mQ,i 2 ,Am 1( /2,m^ 2 and O = (where stands for all the scalar multiplets). The 
dimensionless coefficients Coi and Cqj depend only on the coupling constants at Mqtjt an d on 
t = log 10 (MGUT/^siJSY)- Let us consider the limit of small ji where the correct electro- weak 
symmetry breaking requires that m| becomes negative. In fig. 1 we have plotted the coefficients 
Cqj for O = m\^m\ d and / = m^ 2 ,mg as a function of m t . We see that in the limit rriy 2 3> m\ 
where m\ ~ Cm^mf /2 TO i/2 = ^ ne correc t spontaneous electro-weak symmetry breaking 

can occur already for a top mass of m t <: 60 GeV assuming /i is sufficiently small. We also find that 
m 2 ^ < m\ d due to the different hypercharges of the top and bottom quarks which automatically 
guarantees that tan/? > 1. On the other hand, in the limit m\i 2 <C m 2 the spontaneous electro- 
weak symmetry breaking can only occur for mt ^ 165 GeV. In addition, we see that > m 2 ^ 
which yields tan/3 < 1. Thus, this scenario predicts the wrong hierarchy for the quark masses of 
the third generation and is ruled out. 

Let us now return to the problem of satisfying eq. (5) without fine-tuning. By solving the RG 
equations we obtain 

B ~ C BB B + C B m 1/2 m 1/2 + CbaA , (9) 

where the coefficients are Cbb — 1? C B m 1/2 — ~ 0-8 and Cba — 0.7. Here, we have used cxgut — 
4.2 x 10~ 2 for the unified gauge coupling, ay = y^/An = 3 x 10~ 2 for the top yukawa coupling 
(this corresponds roughly to a top quark mass of mt = 180 GeV; in our numerical work we shall 
instead fix m t = 180 GeV by varying ay) and t = 14. By solving eq. (9) and (7) we obtain Rb 
as a function of A, ro^ 2 and B. 

Note that the definition in eq. (7) is not very satisfying in certain regions of the parameter 
space. In particular, in the region where m 1 / 2 — A we find that Rb — 1 due to a large cancellation 
on the right hand side of eq. (9) (B being negligible). However, this does not mean that no fine- 
tuning is required along this curve in the m^2 — A plane but rather that it has shifted from 
tuning B to tuning A and m^ 2 - Thus, in order to avoid such a cancellation we simply drop A 



from eq. (9) and obtain a more reliable estimate of the amount of fine-tuning 

R B a — ^~ 
B 

As we have argued above, we can derive a lower limit for m l j 2 and thus also for Rg by requiring 
the correct electro- weak symmetry breaking [eq. (4) and (5)]. We obtain the condition 

In Mgut 

<-< = ^2 / {3y 2 X t -3y 2 X b -y 2 X T )dt 

In m z 

= -0.13ml + 0.26m 2 1/2 - 0.04m 1/2 A - 0.01A 2 £ m\ , 

where Xt = m\ + m \ R + m ^ + A 2 etc. Note that the coefficient of m\ in eq. (8) is negative in all 
models where the top and bottom Yukawa couplings are approximately equal due to the additional 
effects of the r Yukawa coupling. This can be compensated by a large gaugino mass parameter 
[16] which enters with a positive sign. However, it requires m^2 = O{. m o) which contradicts an 
approximate R symmetry. 

On the other hand, for the \x parameter we find from eq. (5) for the particular value of ay- 
states above 

fi 2 ~ -\m\ - m 2 hu = -\m\ + 0.06m^ + 2.35m 2 /2 - 0.3m 1/2 A + 0.07 A 2 , (12) 

or ~ 1.5mi/2 which is inconsistent with an approximate PQ symmetry. The situation will not 
change significantly for a different choice of ay- Thus, in this minimal version of the model, where 
all the mass parameters are universal, the large value of tan (5 can not be protected by R or PQ 
symmetry can only be obtained by fine-tuning. 

3. Non-universal mass parameters 

The simplest extension of the minimal SO (10) SUSY-GUT model presented above is the 
inclusion of non-universal mass parameters. This entails a loss of predictability but is certainly 
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favoured over the introduction of additional fields. In SO(10) models gauge invariance puts strong 
constraints on the form of these non-universal terms 

Kon-univ = Am ii F i F i + ( A i3 F i F 3 h + H - C + Amj^/l , (13) 

i,3 i,j 

where Am?. (A?) are hermitian (symmetric) and traceless. [Note, that we have absorbed the 
Yukawa couplings into the definition of Aij.] The terms in eq. (13) are naturally generated by 
integrating out additional fields at Mqut [17; 18] or by RG evolution from Mp to Mgut [H]- 
However, all the terms in eq. (13) will do very little to solve the basic problem of rendering the 
left hand side of eq. (11) positive without a large value of m 1 / 2 - What we need are different mass 
parameters for the two Higgs doublets. In minimal SO(10) symmetric theories, where h u and 
are assumed to be embedded in a single representation, the universality of these two parameters 
is based on gauge invariance and cannot be broken explicitly. Of course, it can be achieved by 
enlarging the Higgs sector. However, as we will point out in the next section such an extension 
is not necessary since the same gauge symmetry that enforces universality of the Higgs mass 
parameter will also break this universality via the U(l)x D-term once the SO(10) symmetry is 
broken. 

3.1 Non-universal masses from D-terms 

Let us assume for simplicity, that the SO(10) gauge group is broken to the SU(3) C <S> SU(2) L 
<8>U(1) Y standard model group at Mgut without any intermediate scale and that the breaking of 
the U(l) x symmetry in SO(10)dSU(5)®U(1) x to SU(5) occurs only via the VEV of #(16) and 
H(16). This can be achieved e.g. with the superpotential 

W = XS(HH - m 2 ) + Ws , (14) 

where m = O(Mqut), S is a SO(10) singlet and Ws contains the quadratic and trilinear term 
in S* In the limit of unbroken R symmetry [i.e. A — B — C — m 1 / 2 = 0) the potential can be 

* Note, that the superpotential in eq. (14) can only be considered as a simple toy model. New fields introduced to 
break SU(5) and to solve the doublet/triplet problem [19;20] can modify our results quantitatively. However, 
the qualitative features of generating a non-vanishing D term will remain. 
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written as V — Vp + Vd + Koft where 



.90 



Koft = 5^ «4 



(15) 



Here the sum is over all SO(10) scalar multiplets = H, H,etc. Universality means that m| = wig 
at Mp for all scalar fields, 0. Note that gauge-invariance requires that all fermions of one generation 
obtain the same mass term at Mqut- With the decomposition H = Hq(10, —1) © ifo(5,3) © 
H N (1, -5) and H = Hq(10, 1)@H d (5, -3)®H n (1, 5) under SO(10)dSU(5)©U(1)x we can write 
the minimum of V in the supersymmetric limit (i.e. V so n = 0) as (Hx) = (-ffjv) — m with all 
other VEVs vanishing. If we include V^ ft then the minimum of V shifts by 0(m^jm) and we will 
find in general that the D-term corresponding to the U(l)x gauge symmetry no longer vanishes 
[21] 

Table 1. The decomposition of F and h under SU(5)<g>U(l)x and SU(3) c <g)SU(2) L ®U(l)y 





F 


h 


SO(10) 


16 


10 


SU(5)©U(1) X 


Q © D © N 


© h D 


SU(5) 


10 


5 


1 


5 


5 


u(i)x 


-1 


3 


-5 


2 


-2 


SU(3) c <g>SU(2) L (g>U(l)y 


q © u c © e c 


d c ® I 


N 


© h u 


*d © &d 


SU(3) C 


3 


3 


1 


3 


1 


1 


3 


1 


3 


1 


SU(2) L 


2 


1 


1 


1 


2 


1 


1 


2 


1 


2 


U(l)y 


l 

3 


4 
3 


2 


2 
3 


-1 





2 
3 


1 


2 
3 


-1 



(D X ) = 9 f (X N (H N y (H N ) + X R (H N y (H N )) = ^ 2 



(m\ - m 2 H ) , 



(16) 



where X R = —Xx = 5 are the U(l)x charges listed in table 1 with the normalization gx = 



Vl/100. The additional non-universal soft SUSY breaking mass terms for the Higgs bosons, 
squarks and sleptons created by (Dx) ^ can be written as 



AV soft = -Am2,^X $ $t$, (17) 

where the sum is now over all SU(5)®U(l)x multiplets listed in table. 1. The non-universal mass 
term is given by [22] 

Am 2 x = g x (D X ) = ^- {m\ - m\) . (18) 

Note that Am x is independent of gx and the parameters of W indicating that this result is quite 
general. 

Of course, in the limit where universality is unbroken the right hand side of eq. (18) vanishes. 
However, the breaking of universality at Mqut due to radiative corrections is almost inevitable. 
For example, let us include interaction terms of the form 

WD±(y H HHti + y s HHh') , (19) 

where h' = h'(10) and h! = h'(10) are different from the Higgs field h. If we impose universality 
(m^, = m 2 H = mg = m^) at Mp and evolve the mass parameters to Mqut we obtain to first 
order [23] 

m 2 H (M G jj T ) = m 2 H (M P ) - ^ (2m 2 H + m 2 ,) In i^A + - , (20) 

>■ GUT / 

where we have omitted all the irrelevant terms. The non-universal mass term then becomes* 

^^ x0 (3). (21) 



★ an interesting alternative to eq. (19) is W D J^^iHFih' . It naturally explains the plus sign of Am|. Here, 
one combination of the three -Dj(5, 3) C i*i(16) (which should lie dominantly in the first two generations) 
combines with the h[j(5, 2) C /i'(10) and becomes massive while the h' D (5, -2) C h'(lO) remains light. This 
might also explain the bad mass predictions for the down-type quarks of the first two generations. 
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Finally, by running the mass parameters to Msusy we find to a good approximation 

A\/ soft ~ -Ami E ( X <h ~ tlfo > ( 22 ) 

i 

where the sum is over all MSSM scalar particles with the charges listed in table 1. In particular, 
we find 

™L-™L-2Am^ + ..., (23) 

where the ellipse stands for the contributions of the universal terms [eq. (11)]. 
3.2 R symmetry 

We now come to a numerical evaluation of the fine-tuning in the minimal SO(10) model 
including non-zero values of Am^. First, we note that for Am^/m^ ~ 7% the right hand side 
of eq. (11) becomes positive independent of the value of m 1 / 2 - By inspecting the RGEs we find 
that the squark mass parameters for the first two generations are numerically irrelevant. The 
same is true for the flavor changing neutral current (FCNC) mediating off-diagonal entries once 
we impose the experimental constraints [24]. Thus, the only important scalar mass parameters 
are sfermion mass parameter of the third generation and the Higgs mass parameter which we take 
to be universal for now (i.e. we set Am| = Am^ = 0). The top-quark mass (or alternatively ay) 
can be fixed by imposing exact r-bottom Yukawa unification. However, there are still significant 
theoretical uncertainties. This becomes manifest if we try to unify the down-type quark and 
lepton Yukawa couplings of the first two generations. Here, one finds that the predicted value 
for (m^m^) / (m s m e ) ~ 1 differs from the experimental value [5] by a factor of order 10. The 
prediction for m^ from Yukawa unification in a minimal model is, therefore, not expected to be 
reliable to more than about the mass of the strange quark m s & 300 MeV [5] which translates 
into a huge incertainty in the prediction of ay. Thus, we allow values of mt to vary within its 
experimental bounds [6] (in all plots where mt is fixed we chose mt = 180 GeV). 

In fig. 2 we present contours of constant fine-tuning, parameterized by Rb = 3, 10, 30, 100, 
300, 1000 in the mj^-Aml plane normalized to mo = 500 GeV for m t = 180 GeV and A = 0. 
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The shaded area is ruled out by requiring that tan (3 > 1. Acceptable values of tan (3 for Am| = 
can only be found for values of m^ 2 — 0.7mo [16] where the fine-tuning is excessive and exceeds 
the lower bound of Rb ^ tan/5 derived in ref. 9 by more than an order of magnitude. However, 
in the region where Am^/m 2 , ~ 7% the value of and thus also the value of Rb can become 
much smaller. 

We now investigate whether the region of acceptable fine-tuning is compatible with the phe- 
nomenological constraints. Here, we only impose the constraints coming from direct particle 
searches in order to keep our results as general as possible. In particular, we do not impose any 
constraints coming from the proton lifetime [25], the decay b — > 57 [26], or relic abundance [27]. 

The most important constraint can be derived from non-observation of the decay of the Z 
boson into a chargino pair (Z — > xfxj, hJ — 1>2) at LEP [5]. It implies that m 1 / 2 ^ 50 GeV 
which is similar to the bound on m^2 from the decay Z — > XiX% hj = 1)2,3,4. This constraint 
alone eliminates almost the entire region in fig. 2 where Rb 30. No significant constraints arise 
from the non-observation of sfermions since this only constrains the region of small m\ (for fixed 
mi 12) where Rb is large anyways. Thus, we expect to be able to decrease -Rb by increasing mo. 
However, in the region where mo m z we have to take into account the fine-tuning required to 
satisfy eq. (4). In the large tan/5 limit we find that \i has to be tuned such that the denominator 
in eq. (4) vanishes as sin 2 2/5. The corresponding fine-tuning is given by 

R _j^_^_ = j^_^_ = 2 j^_ (24) 

M m 2 d\i 2 m\ da 2 m 2 ' 

where we have substituted m? for m 2 . This choice is more appropriate since eq. (4) is derived 
from minimizing the Higgs potential. At tree-level these quantities are equal up to terms of the 
order of 1/ tan 2 /5 but large radiative corrections can increase m^, by a factor of up to two [28]. 
Since we see from eq. (12) n has to grow with uiq (with a coefficient that depends on ay) the 
absence of fine-tuning defines an upper limit of tuq. 

We now proceed to evaluating the combined fine-tuning, parameterized by R = RbR/j, numer- 
ically. Four of the nine GUT input parameters, A, B, /i, m^2, m^, Am^, «GUT) ®y and Mqut, 
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are fixed by experimental observables [5;30] 

sin 2 9 W = 0.231 , 
a£ = 127.9, 

(25) 

m z = 91.187 GeV, 
m T = 1.7771 GeV. 

We fix the non- universal term Am^/mg = 0.1 whose sole purpose it is to assure that m 2 ^ < m 2 ^ 
but which is otherwise of minor importance. Furthermore, we chose m t = 180 GeV, and A = my 2 
and we present the results in the m 1 / 2 ~mQ plane in the ranges 380 GeV & uiq ^ 1 TeV and 
25 GeV & my 2 & 1 TeV. 

Once these parameters are specified we can compute the full supersymmetric particle spectrum. 
We define Mqut as the scale where SU(2)x,<g)U(l)y coupling constants intersect. Here, we include 
all leading logarithmic SUSY threshold corrections by decoupling the superpartners from the 
RGEs at scales below their masses as well as the finite thresholds enhanced by a factor of tan f3 
[31]. The fi parameter is fixed by imposing the minimum conditions [eq. (4) and (5)]. In fig. 3 we 
present contours of (a) constant |//| = 150 GeV, 200 GeV, 300 GeV, 400 GeV, 500 GeV, 600 GeV, 
700 GeV and (b) constant M GUT = 1.75, 2.0, 2.25, 2.5, 2.75, 3.0, 3.25 x 10 16 GeV. We find that 
\jj,\ = 0{my^) unless my^ & m z where we have = 0(mo/3) 

The lightest squarks and sleptons are expected to be the ones of the third generation due to 
the effects of the Yukawa couplings. In fig. 2 we present contours of (a) constant lightest bottom 
squark mass, , (b) constant lightest r slepton mass, Mf i; (c) constant CP-odd Higgs mass, 
m^o, and (d) constant lightest CP-even Higgs mass, m^a in the mo-m^ 2 plane. Here, M^, 
and m^o are evaluated at tree-level and the prediction for m^o includes one-loop and two-loop 
radiative corrections [28;32]. We see that in all of the parameter space under consideration the 
particle spectrum is heavier than the experimental limits from direct particle search at LEP or 
TEVATRON [5]. 

From unification of gauge and Yukawa couplings we can also predict a s and m&. In fig. 5 we 
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present contours of (a) constant a s = 0.126, 0.128, 0.13, 0.132, 0.134 and (b) constant pole mass 
mj = 3.8, 4.0, 4.25, 4.5, 4.75, 5.0, 5.2 GeV in the mj^-mo plane. The values of a s is about one 
standard deviation above the world average value a s (m z ) = 0.12±0.01 [5]. This discrepancy which 
is still within the experimental errors might also be indication for physical effects such as GUT 
threshold corrections. In order to prevent these possible effects to propagate into our prediction 
for mj, we use the experimental value a s (m z ) = 0.112 from deep inelastic scattering [5]*for the 
evolution of mj from m z to m^ and to convert the running mass into the pole mass [33] using 
two-loop threshold corrections and two-loop (3 functions. We find that the so obtained prediction 
of m^ is somewhat smaller than experiment for large values of m^ 2 due to large SUSY threshold 
corrections to m^ [10;18;31]. The cusps in the contours of fig. 5 (a) and (b) for m 1 / 2 — 100 GeV 
correspond to the decoupling of the wino and might be an artifact of our leading logarithmic 
approximation of the threshold corrections. This behaviour is expected to change slightly if one 
uses full one-loop threshold corrections. However, we do not expect any significant effects from 
non-logarithmic terms. 

Finally, in fig. 6 we present contours of constant R = 50, 100, 200, 500, 1000, 2000 in the 
m i/2~ m o plane. We see that the situation has improved significantly with respect to the universal 
case where only the region where m^ 2 ^ 0.7mo i.e. R = O(10 3 ), is allowed. Nonetheless, a lower 
limit of R ^ 50 remains if we impose the experimental constraint ffl^ 2 ^ 50 GeV. 

3.3 PQ symmetry 

In section 3.2 we have investigated the scenario where large values of tan/5 are protected from 
radiative corrections by an approximate R symmetry. However, this symmetry together with 
the experimental lower limits on the chargino masses require a heavy sfermion spectrum [i.e., 
mo = 0(tan(3my2)] which in turn re-introduces a hierarchy problem. The same problem arises if 
we instead try to impose an approximate PQ symmetry. The best option is therefore to impose 
R symmetry and PQ symmetry simultaneously. Such a correlation between R symmetry and 

* we favor this value over LEP data since it is extracted at an energy scale at which the dominant contributions 
to nib arise. 
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PQ symmetry is quite natural in models where the \i parameter is generated dynamically [34]. 
The correct electro-weak symmetry breaking can be achieved in this case by further relaxing the 
universality constraints. Namely, we can fix ji to any arbitrary value by allowing for non-zero 
values for Am^. In this case we have to tune Am^ instead of \x in order to obtain the correct 
value for the VEVs. Thus, the fine-tuning parameter defined in eq. (24) has to be replaced by 



Am^ dm\ ^■ m h 

R H — > R Arnl = 2 ~ 1,2 2~ ' ( 26 ) 



For our numerical analysis we will chose m^2 = — At = 120 GeV which yields a sufficiently 
heavy chargino/neutralino mass spectrum (the sign of \i is chosen such that the SUSY threshold 
corrections to m& are negative). 

Furthermore, in fig 5 (a) we have seen that the prediction for a s tends to be somewhat larger 
than the experimental value. Thus, let us also include the effects of GUT threshold corrections. 
These corrections can be included by modifying the boundary conditions at Mqut 

— OL2 — 013(1 + e) , (27) 

where all the corrections are absorbed in the single parameter e. In fig. 7 we present the prediction 
of a s (dashed curve), m t (dotted curve) and m& (solid curve) as a function of e. We see that for 
e = 1% the prediction of a s , m^ and mt for fixed «gut an d ay decrease by about 3%, 3% and 
1%, respectively. The magnification of the threshold corrections to a s by a factor of a s /aGUT can 
be understood easily. On the other hand, the corrections to and m t are of higher order and 
more complicated. Nonetheless, their dependence on e is quite strong. 

In any realistic model these threshold corrections are indeed expected to be significant [35]. 
Thus, rather than imposing exact unification at Mgut we determine the threshold correction, e 
by fixing a s = 0.12 i.e. we set e = 3%. We can now use the value of a s obtained from unification 
for the evolution of m& also at scales below m z rather than using experimental results as we have 
done in fig. 5 (b). Furthermore, we vary mt within its phenomenologically allowed range. 
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In fig. 8 we present contours of constant Amljm\ in the m$-mt plane. The contour Am| = 
corresponds to the case of universal Higgs and squark masses. We see that in the limit of 
approximate R symmetry and PQ symmetry (in the present case: /i/mo,m^ 2 /mo = 0.4 ~ 0.08) 
and with Am| being the only non-universal soft SUSY breaking parameter the correct electro- 
weak symmetry breaking can only be achieved for mt = 162 ~ 170 GeV. In turn, we can say 
that if the top quark mass falls within this range and the Higgs and top squark mass parameter 
are universal we naturally obtain a large hierarchy m\jm\ <C 1 in the limit of approximate R 
symmetry and PQ symmetry. This point is illustrated in fig. 9 where we present contours of 
constant R in the m t -mo plane and find that for the range of m t denoted above the fine-tuning 
is quite small. Note, the existence of a contour R = in fig. 9 which is equivalent to the contour 
Am^ = in fig. 8 is an unphysical artifact of our definition of the fine-tuning in eq. (26). For small 
values of Am^ it is more appropriate to replace Am^ by the average deviation from universality, 
i.e. 



Am h -M/ 3 



(Ami) 2 + (Ami) 2 ] . (28) 



With this replacement the fine-tuning parameter R would remain roughly constant within the 
contours \Aml/m^\ = Am\jm\ = 1/10. It is interesting to note that these contours correspond 
roughly to \R\ = 50 ~ tan/3, a limit that has already been established by more general consider- 
ations in ref. 9. Note also that R is largely independent of mo but exhibits a strong dependence 
on mt. 

In fig. 10 we present contours of constant pole mass m& in the m^-mo plane. Here, the strong 
coupling constant lies in the range a s = 0.120 ± 0.002. We find that the predicted value of mi 
lies outside the experimental error bars (shaded area). The reason for the difference between this 
plot and fig. 5 (b) is twofold. First we have used a larger value of a s from m z to m&. Second, 
by imposing R symmetry and PQ symmetry we have severely constrained the size of the SUSY 
threshold corrections to m^ which are proportional to [xmgjm\ < 1. In particular, we see that 
without any GUT corrections for the Yukawa couplings there is no overlap between the area 
of small fine-tuning (say \R\ ~ 50) and the area where mj, falls within the experimental error. 
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However, as we have pointed out earlier additional (possibly non-renormalizable) terms have to be 
introduced to explain the masses of the first two generations. If we keep in mind that in minimal 
SO(10) the universal Yukawa coupling determined by the charm quark mass m c = 1.2 GeV can 
only account for a small fraction of the strange quark mass then we can expect the theoretical 
uncertainty in m;, to be of the order of m s & 300 MeV [5]. Thus, by allowing the entire region 
where the predicted values of in fig. 10 is below 5.2 GeV there will be a region in parameter 
space with m t ~ 165 GeV where the fine-tuning is of the order of tan (3. 

4. Conclusions 

We have presented a detailed analysis of the naturalness of the large tan (5 region required in 
minimal SO (10) SUSY-GUT models. 

We confirm that in the case of universal soft SUSY breaking parameters the particle spec- 
trum is dominated by a large gluino mass and the fine-tuning needed for the correct electroweak 
symmetry breaking is (9(10 3 ) and excessive. 

The situation improves if we include a non-zero U(l)x -D-term. This opens a whole new 
range in parameter space where R symmetry is approximately conserved and the required fine- 
tuning is only (9(tan,3). The region of smallest fine-tuning is characterized by a light chargino 
and neutralino just above the experimental bounds. 

Finally, we use non-universal soft SUSY breaking parameters for the Higgs sector and the 
sfermion sector in order to impose an additional approximate PQ symmetry. We find that in this 
scenario the top quark mass is the crucial parameter. In particular, for mt = 0(165 GeV) almost 
universal boundary conditions will result in a hierarchy between the electroweak scale, m z , and 
the SUSY scale, mo. This can further improve the fine-tuning situation in the large tan/? region. 

This scenario is characterized by a m t significantly below the its IR fixed point. As a result, the 
prediction of m^ is above the experimental value. However, we argue that the resulting prediction 
falls still within the significant theoretical errors related to GUT scale threshold correction and to 
the origin of the Yukawa couplings for the first two generations. 
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FIGURE CAPTIONS 

1) The coefficients Coi (O = m\^m\ d and I = mQ,rriy 2 ) as a function of m t 

2) contours of constant Rb = 3, 10, 30, 100, 300, 1000 in the mi^-Am| plane (normalized to 
m ) 

3) contours of (a) constant = 150, 200, 300, 400, 500, 600, 700 GeV and (b) constant 
M GUT = 1.5, 1.75, 2.0, 2.25, 2.5, 2.75, 3.0, 3.25 x 10 16 GeV 

4) contours of (a) constant M~ bi = 150, 200, 300, 400, 500, 750 GeV, 1 TeV, (b) constant 
M fl = 250, 300, 400, 500, 600, 750 GeV, (c) constant m A o = 60, 100, 150, 200, 300 GeV 
and (d) constant to^o = 105, 110, 115, 120, 125, 130, 135 GeV in the mo-m^ plane. We 
set A = my 2 and m t = 180 GeV 

5) contours of (a) constant a s = 0.126, 0.128, 0.13, 0.132, 0.134 and (b) constant pole mass 
m b = 3.8, 4.0, 4.25, 4.5, 4.75, 5.0, 5.2 GeV in the m 1/2 -m plane for m t = 180 GeV and 
A = m 1 / 2 

6) contours of constant R = 50, 100, 200, 500, 1000, 2000 for m t = 180 GeV and A = m 1/2 

7) m t , mi, and a s as a function of the GUT threshold corrections parameterized by e 

8) contours of constant Abi^/toq in the m§-m t plane. We fix = ~H = 120 GeV for an 
approximate R symmetry and PQ symmetry 

9) contours of constant \R\ = 0, 50, 100, 200, 500, 1000, 2000 in the m 1/2 -m plane for 
m t = 180 GeV and A = 

10) contours of constant pole mass m& = 4.75, 5.0, 5.25, 5.5 GeV in the m,Q-mt plane 
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Fig. 2. contours of constant i?e = 3, 10, 30, 100, 300, 1000 in the mi/2-Am| plane 
(normalized to mo). 




Fig. 3. contours of (a) constant \n\ = 150, 200, 300, 400, 500, 600, 700 GeV and (b) 
constant M GUT = 1.5, 1.75, 2.0, 2.25, 2.5, 2.75, 3.0, 3.25 x 10 16 GeV 
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Fig. 4. contours of (a) constant M~ bi = 150, 200, 300, 400, 500, 750 GeV, 1 TeV, (b) 
constant M fl = 250, 300, 400, 500, 600, 750 GeV, (c) constant m A o = 60, 100, 150, 200, 
300 GeV and (d) constant m h o = 105, 110, 115, 120, 125, 130, 135 GeV in the m -m 1/2 
plane. We set A = m^j and m t = 180 GeV 
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Fig. 5. contours of (a) constant a s = 0.126, 0.128, 0.13, 0.132, 0.134 and (b) constant 
pole mass mj = 3.8, 4.0, 4.25, 4.5, 4.75, 5.0, 5.2 GeV in the mi^-mo plane for ra t = 
180 GeV and A = m 1/2 
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Fig. 6. contours of constant R = 50, 100, 200, 500, 1000, 2000 for m t = 180 GeV and 

A = m 1 / 2 




Fig. 7. m t , nib and a s as a function of the GUT threshold corrections parameterized 
by e 
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Fig. 8. contours of constant Anij/mp in the m§-m t plane. We fix m^j = — 
f 20 GeV for an approximate R symmetry and PQ symmetry 
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Fig. 9. contours of constant \R\ = 0, 50, 100, 200, 500, 1000, 2000 in the m 1/2 -m 
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Fig. 10. contours of constant pole mass mj = 4.75, 5.0, 5.25, 5.5 GeV in the rao-m t 
plane 



